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Abstract 



In this paper we consider the role that numerical computations - in par- 
ticular Galerkin approximations - can play in problems modelled by the 
3d Navier-Stokes equations, for which no rigorous proof of the existence of 
unique solutions is currently available. We prove a robustness theorem for 
strong solutions, from which we derive an a posteriori check that can be ap- 
plied to a numerical solution to guarantee the existence of a strong solution 
of the corresponding exact problem. 

We then consider Galerkin approximations, and show that if a strong 
solution exists the Galerkin approximations will converge to it; thus if one 
is prepared to assume that the Navier-Stokes equations are regular one can 
justify this particular numerical method rigorously. 

Combining these two results we show that if a strong solution of the exact 
problem exists then this can be verified numerically using an algorithm that 
can be guaranteed to terminate in a finite time. 

We thus introduce the possibility of rigorous computations of the solu- 
tions of the 3d Navier-Stokes equations (despite the lack of rigorous existence 
and uniqueness results), and demonstrate that numerical investigation can 
be used to rule out the occurrence of possible singularities in particular ex- 
amples. 
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1 Introduction 



The Navier-Stokes equations are the fundamental model of fluid flow. How- 
ever, no proof of the long-time existence of unique solutions of these equa- 
tions is currently available. The importance of this long-standing problem 
has recently been highlighted by its inclusion as one of the Clay Foundation's 
Millennium Problems. 

Despite this mathematical obstacle, the equations are used as the basis of 
numerical calculations both in theoretical investigations of turbulence and for 
industrial applications. This paper addresses the validity of such numerical 
computations given the lack of appropriate rigorous existence and uniqueness 
results. 

One might expect that without the existence of a unique solution there 
is no hope of guaranteeing that a numerical approximation is really an 'ap- 
proximation' in any meaningful sense, since it is not clear what is being 
approximated. However, here we prove three results that enable numerical 
and exact solutions to be related in a rigorous way. 

First we provide an explicit check that can be applied to a numerical 
solution to guarantee that it is approximating a solution of the exact problem. 
In particular, this check implies the existence of a unique solution of the exact 
problem over the same time interval as the calculation. Thus the validity 
of a numerical solution can be verified rigorously by a simple a posteriori 
condition. 

We then turn to the particular example of Galerkin approximations. We 
show that given the assumption of the existence of a sufficiently smooth 
unique solution of the exact problem, the solutions obtained by the Galerkin 
method will converge to this exact solution: If one is prepared to take on 
trust that the Navier-Stokes equations are a meaningful model of fluid flow, 
numerical experiments (at least those using the Galerkin method) can be 
justified mathematically. 

Finally, these two results can be combined to show that the existence of 
a sufficiently smooth solution can be verified numerically using an algorithm 
that can be guaranteed to terminate within a finite time (provided that the 
solution exists). 
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Although we concentrate here on spatial discretization (mainly via the 
spectral Galerkin method), it is relatively straightforward to consider fully 
discrete methods, and these are also addressed where relevant. 

In order to discuss our results in a little more detail, we introduce here 
the model and its abstract setting. We will consider the 3d Navier-Stokes 
equations for an incompressible fluid 

du ~ 

— - vku + (u ■ V)m + Vp = f(t) with V ■ u = (1) 
at 

on a periodic domain Q = [0, L] 3 with the additional (convenient) restriction 
that fgii = Jq f = 0. However, we would expect similar results to hold 
for more physically realistic boundary conditions (e.g. Dirichlet boundary 
conditions on a bounded domain). 

In order to recast the equations in their functional form (for full details see 
the monographs by Temam (1977), Constantin & Foias (1988), or Robinson 
(2001)) we let J4? be the collection of all divergence- free smooth periodic 
vector-valued functions on Q with zero average, and set 

H = closure of in [L 2 (Q)} 3 
V = closure of in [H^Q)} 3 . 

We use | • | and (•, •) for, respectively, the norm and inner product in [L 2 (Q)} 3 . 

Denoting by II the orthogonal projection of [L 2 (Q)] 3 onto H, we apply IT 
to (Q) and obtain 

— + vAu + B(u,u)=f(t), (2) 

where A is the Stokes operator Au = —II Am (in fact Au = —Am when 
u G D(A) for the periodic case which we consider here), B(u, u) = H[(u- V)w], 
and / = II/. The pressure term disappears since gradients are orthogonal 
(in L 2 ) to divergence-free fields. 

We note that in the space periodic case we have D(A m ^ 2 ) = H m D V; for 
simplicity we denote this space by V m , and equip it with the natural norm 

||n|| m = \A m ' 2 u\ 

which is equivalent to the norm in the Sobolev space H m . 
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In section 2 we present a very simple ODE lemma that forms the basis of 
what follows. Section 3 recalls some classical estimates on the nonlinear term 
B(u,u) in the spaces V m (m > 3), and reproduces the proof of a classical 
regularity result. 

Section 4.1 proves the robustness of strong solutions: We show in the- 
orem El for m > 3 that if the initial data uq G V m and forcing f(t) G 
L 2 (0, T; V 171 ' 1 ) fl L 1 (0, T;V m ) give rise to a strong solution, then so do 
'nearby' initial data Vq G V m and forcing functions g G L 2 (0, T; V" 1 ^ 1 ), 
under the explicit condition that 



T 

\"(> - ''nil,,, + / - g(s) 

10 



ds 



< 



exp 



r 







\u(s)\\ m+1 ) ds 



(note that this depends only on the solution u). A similar result holds for 
the Euler equations under slightly modified hypotheses, and this is shown in 
section 4.2. 

As corollary 03 in section 5 we use the observation that any sufficiently 
smooth numerically computed solution u is the exact solution of the Navier- 
Stokes equation for some appropriate forcing function to turn the previous 
robustness result into an a posteriori test for the existence of a strong so- 
lution based on numerical observations. Leaving the precise definition of 
'sufficiently smooth' to the formal statement of the result, we show that if u 
is a 'good' approximation to the solution of 



dv 

— + vAv 
dt 

in the sense that 



B(v,v) 



f(t) 



with 



v (0) =v E V n 



(3) 



\u(0) - v \\ m + 



< 



c T 



exp 



du 
~dt { 



c r . 



+ uAu(s) + B(u{s),u(s)) - f(s) 



ds 



U(S 



u(s 



m+1. 



ds 



then v(t) must be a strong solution of © on [0, T] with v G L°°(0, T; V m ) fl 
L 2 (0,T; V m+1 ). Crucially, this condition depends only on the numerical so- 
lution u (and the given initial data and forcing, Vq and f(t)). 
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In section 6 we consider the convergence of approximations to solutions of 
the Navier-Stokes obtained via a Fourier Galerkin method. Here we assume 
the existence of a sufficiently regular exact solution (in fact we take u$ £ V m 
and / £ L 2 (0,T;V m ) and assume the existence of a strong solution), and 
demonstrate in theorem |U] that under this condition, its Galerkin approxima- 
tions u n do indeed converge to the correct limit in both L°°(0, T;V m ) and 
L 2 (0,T; V m+1 ). 

Finally we combine our a posteriori test and the convergence of Galerkin 
approximations in theorem |H] to show that if, for a given choice of initial 
data M e V m and forcing / £ L^O, T; V m ) n L 2 (0, T; \Z m_1 ) (where m > 3), 
a strong solution does exist for some time interval [0, T] , then this can be 
verified computationally in a finite number of steps. (Essentially we show 
that for a sufficiently large Galerkin calculation we can ensure that the a 
posteriori test of corollary El must be satisfied.) 

It should be noted that we do not aim here to prove results that are 
optimal with regard to the regularity of solutions, but rather to take sufficient 
regularity (initial conditions in the Sobolev space H m with m > 3) that 
the arguments are at their most straightforward. Similar results for initial 
conditions in H 1 and H 2 are possible, and will be presented elsewhere. 



2 A simple ODE lemma 

The following simple lemma, after that in Constantin (1986), is central to all 
the results that follow. Although all solutions of the equation 

y = S(t) + ay 2 with y(0) = yo > and a > 

blow up in a finite time, if yo and 5(t) are sufficiently small (in appropriate 
senses) then the solution can be guaranteed to exist on [0, T]. 

Lemma 1. Let T > and a > be constants, and let S(t) be a non- 
negative continuous function on [0,T]. Suppose thaty satisfies the differential 
inequality 

%<m + uy 2 with y(0)=y >0 (4) 
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and define 

V = Vo+ / S(s) ds. 
Jo 

Then 

y(t) < y 3 — (5) 

1 — arjt 

while < t < T and the right-hand side is finite. In particular if 

ar]T < 1 (6) 

then y(t) remains bounded on [0, T], and clearly y{t) — > uniformly on [0, T] 
as r] — ► 0. 

Proof. Observe that y(t) is bounded by Y(t) where 
dY 



aY 2 with Y(0)=y + 5(s)ds, 
dt J 

which yields (0) and the remainder of the lemma follows immediately. □ 



3 A classical regularity result 

In this section we prove a slight variant of a classical regularity result which 
can be found as theorem 10.6 in Constantin & Foias, 1988. 

First, we recall the following bounds on the nonlinear term, which we will 
use repeatedly (see Constantin & Foias (1988) for a proof): 

\\B(u,v)\\ m < c m \\u\\ m \\v\\ m+1 m > 2, (7) 

\(B(w,v),A m w)\<c m \\v\\m+i\\w\\l m>2, (8) 

\(B(v } w) } A m w)\<c m \\v\\ m \\wf m m>3. (9) 

Theorem 2. Let u G L°°(0,T; V) H L 2 (0,T; V 2 ) be a strong solution of the 
3d Navier-Stokes equations with uq e V m and f G L 2 (0, T; V" 1 ' 1 ). Then in 
fact 

u G L°°(0, T; y m ) n L 2 (0, T; V m+1 ) 

and 

du/dteL 2 {{],T-V m - 1 ). (10) 
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Proof. We give a formal argument which can be made rigorous using the 
Galerkin procedure. The proof is inductive, supposing initially that u G 
L 2 (0, T; V k ) for some k < m. Taking the inner product of the equation with 
A k u we obtain 

^lkll' + Hhl| 2 fe+ i<l(^K«)^Ml + (/,^) 

and so using (|SJ) (valid here for k > 2) 
1 d 



2dt" M '^ + Z/ " M "^+ 1 - Cfc IHIfclMU+i + ll/IU-ill M IU+i- 



Therefore 



f t \\u\\i + ^\\u\\i +1 <^\\u\\t + ^^, (11) 



Dropping the term j we have 

d ,, ll9 . / ct\\u 



2 I ' llfc-1 



It now follows from the Gronwall inequality that our assumption u G £ 2 (0, T; V ) 
implies that u G L°°(0, T; V k ). 

Returning to (jllj) and integrating between and T we obtain 

"/ IK«)lll + i^<IKo)||l+^ / \Hs)\\ids+- [ ||/( 8 )||J_id*, (12) 

Jo v Jo v Jo 

which shows in turn that u G L 2 (0,T; V k+1 ). 

Since by assumption u G L 2 (0,T; I/ 2 ), the first use of the induction re- 
quires k = 2, for which inequality (jHJ is valid: we can therefore conclude by 
induction that u G L°°(0,T; V m ) n £ 2 (0,T; \/ m+1 ). Finally, since 

dw 

— — = —vAu — B(u, u) + f 
dt 

and each term on the right-hand side is contained in L 2 (0, T; V™ 1-1 ), the 
bound on the time derivative in (jlOj) follows. □ 



We note here that it follows from this theorem that if Uq and / are smooth 
(in V m for all m) then so is the solution. Different techniques (due to Foias 
& Temam, 1989) can be used to show that the solution is analytic in the 
space variable (i.e. in a certain Gevrey class) provided that the data is. 



8 



4 Robustness of strong solutions 



4.1 The Navier-Stokes equations 



Using lemma ^ we show that if the 3d Navier-Stokes equations have a suf- 
ficiently smooth strong solution for given initial data uq and forcing / then 
they also have a strong solution for close enough data. The argument is based 
closely on that in Constantin (1986) which, given the existence of a strong 
solution of the Euler equations, deduces the existence of strong solutions for 
the Navier-Stokes equations for small enough v > 0. (The same argument is 
used in Chapter 11 of the monograph by Constantin & Foias (1988)). 

We choose to state our primary result for sufficiently smooth strong solu- 
tions, namely those corresponding to initial data in V m with m > 3. This en- 
ables us to use all the inequalities and thereby obtain a relatively simple 
'closeness' condition in (|13p. A similar approach works with strong solutions 
that have the minimal required regularity (u G L°°(0, T; V) fl L 2 (0,T; V" 2 )) 
but the results are less elegant; these results will be presented in a future 
paper. 

Note that while here we concentrate on the robustness of solutions defined 
on finite time intervals, a result valid for all t > given a solution u for 
which J °° ||it(s)||i ds < oo (i.e. which decays appropriately as t — > oo) has 
been obtained by Ponce et al. (1993) for the particular case / = 0. 

Theorem 3. Let m > 3 and let u be a strong solution of the 3d Navier-Stokes 
equations 

du 

— + vAu + B(u, u) = f(t) with u(0) =u e V m 
dt 

and f e L 2 (0,T;V rn - 1 ) H L\0,T;V m ). Then if g G L 2 (0, T; I/" 1 " 1 ) and 

\\u -Vo\\ m + / \\f(s) - g(s)\\ m ds 
Jo 



< 



1 



c T 



exp 



Cm I ( || ^( s ) || m ~\~ ||w(s)|| m +i ) ds 

o 



(13) 



the solution v of 
dv 



— + vAv + B(v, v) = g{t) with v(0) = v G V r ' 
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is a strong solution on [0, T] and is as regular as u. 



We remark here that throughout this paper we consider only Leray-Hopf 
weak solutions, i.e. weak solutions satisfying the energy inequality starting 
from almost every time. 

Proof. Standard existence results guarantee that v is a strong solution on 
some time interval [0, T*). If T* is maximal then ||w(i)||i — > oo as t —>■ T*; 
clearly we also have ||u(t)|| m — ► oo as t — > T*. We suppose that T* < T and 
deduce a contradiction. 

While v remains strong, the assumption that v 6 V m and that g G 
L 2 (0, T; y m_1 ) allows one to use the regularity results of theoremElto deduce 
that 

v E L°°(0, T'\ V m ) n L 2 (0, T'- V m+l ) 

for any T' < T*; we also have dv/dt e L 2 (0, T'; l/ m_1 ) (cf. the argument in 
the proof of theorem^. It also follows from theorem |2 that the solution u 
enjoys similar regularity on [0, T]. 

The difference w = u — v satisfies 

^ + vAw + B(u, w) + B(w, u) + B(w, w) = f - g (14) 

Lit- 

with w(0) = wo = u — Vq. On [0,T*) we know that w is sufficiently regular 
that 

We can therefore take the inner product with A m w and obtain, using (JHJ) and 

^TllMlm < CmlMLHHIm + c m II u \\ m+1 1| W || ^ + C TO ||w||^ + ||/ - #|| m |ML- 

2 (1/ 

(16) 
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Dividing by ||iw|| m yields 1 



^||^||m — \\f ^ || m ^ C m || M || m || W || m C m || M || || W || m C m || W || m 



< 11/ - g\\m + Cm(\\u\\ m + \\u\\ rn+1 )\\w\\ m + \\w 



1 2 



We multiply by exp(— c m J" * ||w(s)|| m + \\u(s) \\ m +i ds) and consider 



y(t) = ||w(t)|| m exp 



for which we obtain the inequality 



dy 



«0) L + «0 



-i) ds 



where 



< 11/ - g\\m + ay 2 with y(0) = \\u Q -v \ 



a = c m exp 



|u(s)||m + IKs)|| m+ l) ds 



(17) 



In this case the condition (jH)) from lemma ^ becomes JT3J). If this is 
satisfied then y(t) is uniformly bounded on [0, T*); the solution can therefore 
be extended as a strong solution beyond t — T*, contradicting the fact that 
T* < T. It follows that v(t) is a strong solution on [0, T] with the same 
regularity as n. □ 



4.2 The Euler equations {y = 0) 



We note here that the dissipative term vAu plays no direct role in the proof 
of theorem El however, it does enter indirectly via the regularity results of 
theorem El that are required to justify the equality in (fTHjl . It is possible to 
circumvent this via an appropriate mollification (cf. Constantin, E, & Titi, 
1994) and obtain the above result (and those that follow) for solutions of the 
Euler equations {y = 0) as well as those of the Navier-Stokes equations. 

1 Of course, one has to worry here whether ||iw(t)|| m is zero. However, if ||u>(io)||m = 
for some to then w(to) = 0, and then the uniqueness of strong solutions in class of Leray- 
Hopf weak solutions implies that u(t) = v(t) for all t > to, i.e. w(t) = for all t € [to^], 
a contradiction. 
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Theorem 4. Let m > 3 and let u £ L 2 (0, T; V Am+3 ) 6e a solution of the 3d 
Euler equations 

(ill 

— + B(u, u) = f(t) with u(0) = u £ V m+3 

LLC 

where f £ L 2 (0, T; I/™- 1 ) nL^T; V m ). If g £ L 2 (0, T; V™' 1 ) and 

T 



u -v \\ m + / ||/(s) -0(s)|| m ds 



1 

Cm-' 



r 



c m / (||w(s)|| m + \\u(s)\\ 

m+1 







(18) 



then the equation 



— + B(v, v) = g(t) with v(0) =v e V m+ \ 
has a solution v £ L°°(0,T; V m ). 

Proof. Suppose that ||u(t)|| m < oo for alH £ [0, T*) but that ||v(t)|| m — > oo 
as t — > T*. As in the proof of Theorem El we assume that T* < T and 
obtain a contradiction. First, observe that taking T" < T* it follows from the 
regularity theory of Beale, Kato, & Majda (1984) that v £ L 2 (0, T"; 1/ m+3 ). 

Now denote by </? £ C^°(1R 3 ) a standard mollifier supported in the unit 
ball, set v ?e (3 ; ) — £~ 3( p{x/t), and write w e for the mollification of u by convo- 
lution with ip € : 

u e (x) = / <p c (y)u(x -y)dy. 



(Note that the convolution is taken over all of M 3 with u extended periodically, 
but the norms and inner products are still taken over Q.) 

We return to (fT4*)l with v set to zero, 

^ + B(u, w) + B(w, u) + B(w, w) = f - g, 

mollify the equation, 
dw 



dt 



+ b(u, wy + b(w, u y + b(w, wy = (f- g y 



12 



and take the inner product with w e in V m . In this way we obtain with 
all quantities replaced by their mollified counterparts, but with the addition 
of three error terms, 

Vuww 3wuw Owwwi 

where 

Suvw = (B(u, v)\ w £ ) m - (B(u e , v £ ), w e ) m . 

Now, note that we have the pointwise identity (cf. Constantin, E, & Titi, 
1994) 

[(u . V)v} e = (u e ■ V)v e + r e (u, v) - [((« - u e ) • V)(v - v')] 

where 

r e (u, v)(x) — J (p e (y) [(u y (x)- V)v y (x)] dy with u y {x) = u{x—y)—u{x). 
It follows that 

||[(« ■ V)vf - (u* ■ VKIU < \\r e (u,v)\\ m + \\((u - u e ) ■ V)(v - v*)\\ m . 
Since in 3d \u(x) — u(y)\ < C\\u\\h 2 \x — y] 1 ^ 2 , it follows that 
\\ u — ue \\w n < ce 1 / 2 ||n||^m+2 and ||r e (w, v)||frm < ^c 2 ||u||#m+2||t>||^m+3j e, 
and so we have 

\8uvw\ < ^||^ e ||,f/ m + 2 ||^ e ||// m + 3 ||w e ||i/m^e. 
Now, use the fact that ||m £ ||h s < ll w lli? s f° r e < -^/2, we obtain 

^Iklm < C m ||M|| m ||w e || m + C m ||w|| m+ i||w e || m + C m ||w e ||^ + ||/ - g\\ m 

+ /c / ^||M|| m+ 2||w e || m +3 + \\w e \\ m+2\\ u \\m+3 \\ w || 1| ^ ||m+3) e 

for each < e < L/2. 
Now define (cf. (fT7jl ) 

a 4 = c m exp (c m [ \\u(s)\\ m + \\u(s)\\ m+l 

ds 
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o 



to obtain 



w e (t)\\ m < 



VT',e 



1 — Ot'TJt'^T' 



for all te[0,T') 



provided that 



u -v \\ m + \\f(s) - g(s)\\ m ds 



< 



c T 



7 exp 



|M||m+2||w e ||m+3 + || ^1 m+2 \\ U || m+3 + || 1| m+ 2 || || m+3 



|w(s)||m + ||u(s)||m+l) d ^ 



(cf. Since both u and t> are regular solutions on [0, T'], it follows that 

r-T' 

||w|| m+2 ||w|| m+ 3 + ||w|| m+2 ||w||m+3 + ||w|| m+ 2||w|| m +3ds 



is finite. Given that ||w e || s < || w || s for all e < L/2 we can therefore let e — > 
and obtain the bound 

r) T > 



w (t)\\ m < 



1 — ax'TjT'T' 



for all te[0,T'} 



(19) 



where now 



provided that 



Vt> = \\u -v \\ m + / ||/(s)-</(s) 
'o 



?7t' < — = exp 

On-' 



k(s)||m + IKs)llm+l) ds 



(20) 



Since (fTHj) holds, it follows that ((20)) is verified for all T" sufficiently close 
to T. We can therefore deduce from (|T5|) that 

77t 



k(t)IU < 



for all te[0,T), 



1 — aTflrT 

and we have obtained a contradiction. It follows that v is a regular solution 
on[0,T]. □ 



Similar techniques should be applicable to extend the other results of this 
paper to the Euler case, but from now on we treat only the Navier-Stokes 
equations. 
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5 Deducing the existence of a strong solution 
via numerics 



As an application of theorem El- and this may be the most significant result 
in this paper - we give an a posteriori test to determine whether or not a 
numerical solution of the 3d Navier-Stokes equations is meaningful. That is, 
we give a criterion depending only on the numerically computed solution that, 
if satisfied, guarantees that the exact equation being approximated possesses 
a strong (and hence unique) solution. Note that in the statement of the 
theorem the provenance of the function u is irrelevant. 

Corollary 5. Let f G L 2 (0, T; V m ) and u G C°([0, T]; V m ) n L 2 (0, T; V m+1 ) 
with 

— + uAu + B(u,u) G l}{U,T-V m )r\L 2 {U,T-V r ' 
at 

for some m > 3. If v £ V m and 



rm—l~- 



KO) - V \\ m + 



^( s )+uAu(s) + B(u(s),u(s))-f(s) 
at 



ds 



< 



cxp 



C r> 



21) 



(IKs)|| m + IKs)|| m+1 ) d S 

then the solution of the Navier-Stokes equation 

fit) with v(0) = v G V m (22) 
is a strong solution on [0,7] with v G L°°(0,T; V m ) n £ 2 (0,T; V m+1 ). 



— + vAv + B(v,v) 
dt 



Proof. The function u is regular enough that it is the (unique) strong solution 
of the Navier-Stokes equation (for the unknown u) 

du du 

— — h vAu + B(u, u) = — — h vAu + B(u, u) with u(0) = u ; 
dt ^dt / 

notional forcing 

note that the conditions on u ensure that the right-hand side is an element 
of L 2 (0, T; V" 1-1 ). We now use theoremElto compare u with the solution of 
((221): in this case the condition (fTB*!) to guarantee that v is a strong solution 
is precisely (|2*Tj) . □ 
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If u comes from a discrete time-stepping algorithm, so that its approxi- 
mate values only specified at times t n (with t n+ \ > t n ) then we can 
define a continuous function u via linear interpolation, 

= , tn+1 _ t u n + — — ^— u n +i for t G (t n , t n+1 \. (23) 

"n+1 "n+1 

Eu n G \^ m+2 then certainly m G C°([0, T]; ^ m+2 ) and du/dt G L 2 (0,T; \/ m + 2 ). 
One can now apply the test of corollary |5] to this function u; if the test is 
satisfied this again proves the existence of a strong solution for (}2*2*|) . 



6 Convergence of the Galerkin approxima- 
tions 

We now turn our attention to one particular form of numerical solution, 
namely the Galerkin approximation. We show that given the existence of a 
suitably smooth strong solution, this numerical method is meaningful in that 
the Galerkin approximations can be guaranteed to converge to the strong so- 
lution. Similar results - convergence given the assumption that a strong solu- 
tion exists - are given for finite element methods by Heywood & Rannacher 
(1982), for a Fourier collocation method by E (1993), and for a nonlinear 
Galerkin method by Devulder, Marion, & Titi (1993). 

In some sense the result of Heywood (1982) that the Galerkin approxima- 
tions of a stable solution of the equations converge uniformly on the whole 
time interval is in a similar spirit: properties of the Galerkin method are 
deduced from an assumption on the full equation. Of course, in this context 
it is perhaps more natural to seek conditions under which one can guarantee 
the existence of such a solution given properties of the Galerkin approxima- 
tions, see for example Constantin, Foias, & Temam (1984) and Titi (1987). 
For similar results for time-periodic solutions see Titi (1991). 

We should emphasise again that, in contrast to some related analyses of 
the Galerkin method (e.g. Rautmann, 1980) no assumption is made on the 
regularity of the Galerkin approximations themselves. 
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6.1 The Galerkin approximation 



We denote by P n the orthogonal projection in H onto the first n eigenfunc- 
tions of the Stokes operator, and by Q n its orthogonal complement. Denot- 
ing these eigenf unctions by {wj}'?L 1 , and their corresponding eigenvalues by 
< Ai < A2 < ■ ■ • , we have 

n 00 

P n U = 2_s( U i W j) W j an d QnU = (u,Wj)lUj. 

j=l j=n+l 

Note that if u G V m then 

00 00 

\\QnU\\ 2 m = ^T\(^^)\ 2 <J2 X T\( U ^)\ 2 =Ml 
j=n+l j=l 

and clearly Q n u — > in V m as n — » 00. 

The Galerkin approximation of (J2J) is obtained by projecting all terms 
onto the space P n H: 

dti 

-^ + vAu n + P n B(u n ,u n ) = P n f(t) with u n (0) = P n u . (24) 



6.2 Convergence of the Galerkin approximation 

Again we present our result for sufficiently strong solutions. With some care 
one can combine the approach of Devulder et al. (1993) with that used here 
to give a proof for strong solutions with minimal regularity; this will be 
presented elsewhere. 

Theorem 6. Let u G V m with m>3,fe L 2 (0,T; V m ), and let u{t) be a 
strong solution of the Navier-Stokes equations 

dw 

— + vAu + B(u, u) = fit) with u(0) = u . (25) 

Denote by u n the solution of the Galerkin approximation \24\j - Then u n — > u 
strongly in both L°°(0,T; V m ) and L 2 (0,T; V m+l ) as n -> 00. 
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Proof. The key, as with the robustness theorem, is to arrange B(u,u) - 
B{u n , u n ) so that it only involves u n in the form u n — u. Writing w n = u — u r 
yields the equation (cf. JUJ) 



dw n 
~dt 



vAw n + P n B(u, w n ) + P n B(w n , u) + P n B(w n , w n ) = Q n f - Q n B(u, u). 



Taking the inner product of this equation with A m w n we obtain 



1_|| ||2 , II u2 ^ || || || II 2 I || || || ||2 _|_ || I ■ 3 

^ dt n n Cm || M || m || if n || "T C m || It || m +l || t^n \\ m ~r C m || 10 n || m 



+ \\Qnf ~ QnB(u,u)\\ m \\w r , 



(26) 



Dropping the term Hl^nJIm+i an d dividing by ||u> n ||m yields 



^||w n || m < c m (||u|| m + ||m|U+i)||w„|| + c m ||w„||^ + \\Q n f - Q n B{u,u) 



Setting 



y n (t) = ||w n (t)|| m exp 



-''»' / (\\u(s)\\ m + \\u(s)\\ m+1 ) ds 





we obtain 

Vn < \\Qnf - Q n B(u,u)\\ m + ay 2 n with 2/ n (0) = ||QnUo| 
where as in the proof of theorem El 



a = c m exp 







Noting that a is independent of n, and that y n (t) is proportional to 
convergence of the Galerkin solutions will follow from convergence 



n \"J || mi 



of y n (t) to zero. Using lemma ^ this will follow from 



\\Q n uo\\m+ \\Qn[f(s) - B(u(s),u(s))]\\ m ds -> as n -> oo, 
./q 



which we now demonstrate. 
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That HQn^ollm as n — > oo is immediate from the definition of Q n . 
To show convergence of the integral term, observe that since 

\\B{u,u)\\ m < c m ||u|| m ||u|| m+1 (27) 

and the regularity result of theorem ^guarantees that u G L 2 (0, T; V m+l ), it 
follows that 

f(s) - B(u(s), u(s)) G V m for a.e. s G [0, T]. 

We therefore know that ||<3n[/(s) — B(u(s), u(s)))\\ m converges pointwise to 
zero for a.e. s G [0, T], while it is clear that 

\\Q n [f(s)-B(u{s),u(s)]\\ m < \\f(s)-B(u(s),u(s))\\ m for a.e. sG [0,T], 

and the right-hand side is an element of L 1 (0, T). It follows from the Lebesgue 
dominated convergence theorem that 

/ \\Q n [f{s)-B{u{s),u{s))]\\ m ds^0 as n -> oo. (28) 

Returning to (|2*Bj) and integrating between and T one can easily show 
the convergence of w n to zero in L 2 (0, T; V m+1 ). □ 

It is clear that if the Galerkin approximations converge then one can 
devise a fully discrete method that converges in a similar sense, and this can 
easily be made precise. In the following theorem one can take u n> s to be the 
linear interpolant of a discrete set of values u nj s(tj) G V m as in . 

Corollary 7. Let u G V m with m>3,fe L 2 {0,T; V m ), and let u{t) be a 
strong solution of the Navier-Stokes equations H3J). For each n, let {u nt s}s>o 
be a collection of functions from [0,T] into P n H with the property that 

u n ,s(t) -> Unit) in L°°(0,T;\/ m ) as 5^0, (29) 

where u n is the solution of the Galerkin approximation \2$ - Then there exist 
5 n > such that 

sup \\u nt s — u\\l°°(o,t-,v™) + \\u n ,s ~ w||L 2 (o,T ; y m + 1 ) — >■ 

0<5<5„ 

as n — > oo . 
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Proof. Since all norms on any finite-dimensional space are equivalent (in 
particular those on P n H = P n V m = P n V m+1 ) the convergence in (J29)) also 
implies convergence in L°°(0,T; V m+1 ), and so in L 2 (0, T; V m+1 ). Otherwise 
the result is immediate. □ 



7 Guaranteed numerical verification of the 
existence of a strong solution 

Our final result combines corollaryEJand theoremElto show that the existence 
of a smooth strong solution can be verified algorithmically in a finite time 
by a sufficiently refined numerical computation. Namely, we show that for 
n sufficiently large the Galerkin solution will satisfy the regularity test (j2Tj) 
from corollary E| 

Theorem 8. Suppose that for some m > 3, Uq G V m , f E L 1 (0,T; V rn ) H 
L 2 (0, T; V" 1 ' 1 ), and that u is a strong solution of the Navier-Stokes equa- 
tions Then there exists an N such that the solution u n of the Galerkin 
approximation \2$ satisfies condition \21\) for every n > N, i.e. will pass 
the a posteriori test for the existence of a strong solution with data (uq, f). 



Proof. First, we note that the convergence of the Galerkin approximations to 
u that is provided by theorem |H1 shows that J Q 
is bounded independently of n. 



\U n (S) 



+ \\uJs) 



m+1 



) ds 



We need, therefore, only show that 



HQn^ollm + 



du r 



dt 



[s) + vAu ri (s) + B(u n (s), u n (s)) - f(s) 



ds 



(the left-hand side of (}2*TjO tends to zero as n — » oo. The requirement on the 
initial condition is trivially satisfied, so we consider here only the 'remainder 
term' m 

du 

""' (s) + vAu n (s) + B(u n (s),u n (s)) - f(s) ds. 



dt 



Now, since u n satisfies the Galerkin approximation (|24p we have 



du 



dt 



-(s) + h>Au n (s) + B(u n (s), u n (s)) - f(s) = Q n [B(u n (s), u n (s)) - f(s)}, 
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and so the remainder term is in fact equal to 



/ \\Q n [B(u n (s),U n (s)) - f( S )]\\ m ds 



(30) 



In particular this shows that solutions of the Galerkin approximation have 
the regularity required to ensure that the integral on the left-hand side of 
(f2Tj) is well-defined. 

The integral in (|30J) coincides with the expression in (|25jl . whose conver- 
gence to zero we showed above in the proof of theorem |Hl except that the 
argument of B is u n rather than u. However, 

B(u n (s), u n (s))-B(u(s),u(s)) = B(u n (s)-u(s),u n (s))+B(u(s),u n (s)-u(s)) 



Since u n -> u strongly in both L°°(0,T; V m ) and L 2 (0,T; V m+1 ) the result 



Once more it is possible to treat fully discrete schemes within a similar 
framework. Supposing as above that a scheme gives rise to a discrete set of 
values u nt s{tj), one can define a linear interpolation u n _s- We then have the 
following result: 



Corollary 9. Suppose that for some m > 3 ; u 6 V m , f G L 1 (0,T; V m ) H 
L 2 (0, T; V m ~ l ), and that u is a strong solution of the Navier-Stokes equations 



and so 




follows. 



□ 
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\29(l. For each n, let {u Ut s}s>o be a collection of functions from [0, T] into 
P n H with the property that 

du n , s (t)/dt^du n (t)/dt in L°°(0,T;V m ) as 5^0, (31) 

where u n is the solution of the Galerkin approximation \2J$ . Then there exists 
an N and a sequence 5 n such that the interpolant u ny s satisfies condition (H7J) 
for every n > N and 5 < 5 n , i. e. the fully discrete numerical solution will pass 
the a posteriori test for the existence of a strong solution with data (uq, f). 

Proof. The interpolant u n ^ satisfies the Galerkin approximation (}2lj) except 
for an error 

_ dUn _ du n: s 
rn ' S ~ dt dt 5 

this error converges to zero in L°°(0,T; V m ) because of (jHlj) (cf. Higham & 
Stuart, 1998). It follows that we now obtain 

T 

\\Q n [B(u n (s),u n (s)) - f(s) +r nj(5 (s)]|| m ds. 

rather than (|30J). where r n j converges to zero in L°°(0, T; V m ). This is clearly 
sufficient to follow the argument in the proof of theorem |HJ □ 



8 Conclusion 



Despite the lack of a guarantee that unique solutions exist for the three- 
dimensional Navier-Stokes equations, we have shown that it is possible to 
perform 'rigorous' numerical experiments. In particular, we have given an a 
posteriori test that, if satisfied by a numerical solution, guarantees that it 
approximates a true strong solution of the Navier-Stokes equations. Remark- 
ably, the existence of such a solution can be verified using such numerical 
computations in a finite time; some computations along these lines will be 
reported in a future publication. 
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